THE GROSS-PITAEVSKII FUNCTIONAL WITH A RANDOM 
BACKGROUND POTENTIAL AND CONDENSATION IN THE 
SINGLE PARTICLE GROUND STATE 

FREDERIC KLOPP, BERND METZGER 

Abstract. For discrete and continuous Gross-Pitaevskii energy functionals 
with a random background potential, we study the Gross-Pitaevskii ground 
state. We characterize a regime of interaction coupling when the Gross- 
Pitaevskii ground state and the ground state of the random background 
Hamiltonian asymptotically coincide. 



1. Introduction 

The purpose of the present paper is to study some aspects of condensation in 
the ground state of the Gross-Pitaevskii energy functional with a disordered 
background potential. As they can be treated very similar, we consider the 
discrete and the continuous setting simultaneously. 

The continuous setting: 

In M.'^, consider the cube A^, = [— L, L]*^ of side length 2L and volume 
|Al| = (2L)'^. In Til ■■= ^^(Al), on the domain Vl := i?^(Ai,), consider 
L ~ + ^^)al continuous self-adjoint Anderson model on A^ 
with periodic boundary conditions. We assume 

• A = Yl'j=i ^"j continuous Laplace operator; 

• 14; is an ergodic random potential i.e. an ergodic random field 
over that satisfies 

VaGN^ ||||a°K,||,,ooL,oc < +00 

where || • ||a;^oo (resp. || • ||i^,oo) denotes the supremum norm in x 
(resp. oj). 

These assumptions are for example satisfied by a continuous Anderson 
model with a smooth compactly supported single site potential i.e. if 



where u G Co°(M'^) and ((x!^)^gA^ are bounded, non negative identically 
distributed random variables. 
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The discrete setting: 

On the finite discrete cube A/, = [—L,LY' r\'L'^ d the cube of side 
length 2L + 1 and volume |Al| = (2L + 1)'^, let H^_^ = (-A + V^)^^ the 
discrete Anderson model on = T-Ll '■= ^^(Al) with periodic boundary 
conditions. We assume 

• (— A)^^ is the discrete Laplacian; 

• 14; is a potential i.e. a diagonal matrix entries of which are 
are given by bounded non negative random variables, say co = 

For the sake of definiteness, we assume that the infimum of the (almost sure) 
spectrum of H^^ be 0. We define 



Definition 1 (Gross-Pitaevskii energy functional [GPEF]). The (one- 
particle) Gross-Pitaevskii energy functional on the cube Al (in the discrete or 
in the continuous) is defined by 



(1) £S,l[^] = KLV,^) + m^\\t 



for if £ Vl and U is a positive coupling constant. 

For applications, it is natural that this coupling constant is related to \A.l\. We 
refer to the discussion following Theorem 3 for details. One proves 

Proposition 2. For any cj € and L > 1, there exists a ground state ip'^^ 
i.e. a vector 99^^ G Vl such that ||v?*^^||2 = 1 minimizing the Gross-Pitaevskii 
energy functional, i.e. 

(2) ^Z = £S'l[^'''']= 

\\<fi\\2 = l 



The ground state (^^^ can be chosen positive; it is unique up to a change of 
phase. denotes the ground state energy of the discrete Gross-Pitaevskii 

functional. 

The proof in the continuous case is given in ['2()]; the proof in the discrete case 
is similar. 

Let and H^j^ respectively denote the Neumann and Dirichlet restrictions 
of to Ax,. Our main assumptions on the random model are: 
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(HO) Decorrelation estimate: the model satisfies a finite range decor- 
relation estimate i.e. there exits R > such that, for any J G N* and 
any sets {Dj)i<j<j, if 

inf dist{Dj,Dj,) > R, 

then the restrictions of V^jj to the domains Dj, i.e. the functions 

)i<i<J) ^I's independent random fields. 
(HI) Wegner estimate: There exists C > such that, for any compact 
interval / and • G {P, N, D}, 

E[tr(lK^:,J)]<C|/|L'^; 

(H2) Minami estimate: There exists C > such that, for I a compact 
interval and • e {P,N,D}, 

f[{H'j^ has at least two eigenvalues in /}] < C{\I\L'^)'^; 

(H3) Lifshitz type estimate near energy 0: There exist constants 
C > c > such that, for L > 1 and any parallelepiped Pl = Ii x ■ ■ ■ x Id 
where the intervals {Ij)i<j<d satisfy L/2 < \Ij\ < 2L, one has 

^g-L<*/c ^ p[{_ff^^^ has at least one eigenvalue in [0, 
P[{i/^p^ has at least one eigenvalue in [0,^^^]}] < Ce"^''/^ 

where Hj^^p^ (resp. H^^p^) is the Dirichlet (resp. Neumann restriction) 
oiH^ to^ Pl. 

Let us now discuss the validity of these assumptions. 

The decorrelation assumption (HO) is satisfied for the discrete Anderson model 
described above if the random variables (w^)^gzd are i.i.d. (HO) clearly allows 
some correlation between the random variables. For the continuous Anderson 
model, it is satisfied if the single site potential has compact support and the 
random variables are i.i.d. 

Under the assumption that the random variables are i.i.d and that their distri- 
bution is regular, it is well known that the Wegner estimate (HI) holds at all en- 
ergies for both the discrete and continuous Anderson model (see e.g. [17, 37, 7]). 
The Minami estimate (H2) is known to hold at all energies under similar reg- 
ularity assumptions for the discrete Anderson model (see e.g. [28, 3, 13, 5]) 
and for the continuous Anderson model in the localization regime under more 
specific assumptions on the single site potential (see e.g. [C]). 
Finally, the Lifshitz tails estimate (H3) is known to hold for both the contin- 
uous and discrete Anderson model under the sole assumption that the i.i.d. 
random variables be non degenerate, non negative and is in their essential 
range (see e.g. [16, 18, 15]). Though the Lifshitz tails estimate is usually not 
stated for parallelepipeds but for cubes, the proof for cubes applies directly to 
parallelepipeds satisfying the condition stated in (H3). 

The main result of the present paper is 
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Theorem 3 (Condensation in the single particle ground state). Assume 
assumptions (H0)-(H3) hold. Denote by ipo the single particle ground state of 
L (chosen to be positive for the sake of definiteness) and by 'p'^^ the Gross- 
Pitaevskii ground state. 
If for L large, one assumes that 



U = U{L) 



L- 



where 



(3) 



(l + (logL)'^-Vd+e)/rf(logL) 

fr'/^ ^fd<3, 
^-i/^log^ ifd = 4, 
C^/"^ ifd>b. 



and € = in the discrete setting, resp. e > arbitrary in the continuous case, 
then, there exists < r]{L) when L +oo such that 



(4) 



GP\ 



1| >r?(L)}] ^ 0. 



The proof of Theorem 3 also yields information on the size of ri{L) and on the 
probability estimated in (4). Note that the assumption (H1)-(H3) can be re- 
laxed at the expense of changing the admissible size for U. 

To appreciate Theorem 3 maybe some comments about the physical background 
of the Gross-Pitaevskii model, its relationship to Bose-Einstein condensation 
and to known results are of interest. Motivated by recent experiments with 
weakly interacting Bose gases in optical lattices (see for example [4]) the fun- 
damental objects of interest are the ground state density and energy, i.e. 



(5) := min ( 



$ecg)L2{Ai) 
11*11=1 



N 

^{-Ai + V{xi)} + v{\xi-Xj\ 

i=l i<i<j<N 



The optical lattice is modeled by the background potential V as shown in 
Figure 1. Assuming a weak interaction limit of the interaction potential v{x,y), 
the continuous A^-particle Gross-Pitaevskii energy functional 

(6) = min / {N\V^{x)\'^ + NV\ip{x)\'^ + 4N'^TTna\^{x)f) dx, 



ipeL^(Ai) J A 

\M\2 = i 

is a mean field approximation of the ground state energy (5), e.g. in three 
dimensions one has 

£-QM 
lim --=rp- = 1. 
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Figure 1. An example of a background potential modeling an 
optical lattice [29] 

( see for example [20] ) The discrete Gross-Pitaevskii model is then a tight 
binding approximation of the continuous one-particle Gross-Pitaevskii func- 
tional [31, 32] 

£^^W]= I (|V(^(x)|2 + y|(^(x)|2 + 4iV7r^a|¥.(x)|^)(ix. 

Another way to derive the discrete Gross-Pitaevskii model starts with a dis- 
cretization of (5) yielding the standard description of optical lattices using the 
Bose-Hubbard-Hamiltonian 

H = - ^ cicn' + "^{crVn - n)nn + ]^U^nl 

\n-n'\=l n 

where Cn, c„, are bosonic creation and annihilation operators and n„ gives the 
particle number at site n (see the survey article [4] and references therein). 
A mean field approximation then yields the discrete Gross-Pitaevskii energy 
functional [2.')]. 

One motivation to study Bose gases is Bose-Einstein condensation, i.e. the 
phenomena that a single particle level has a macroscopic occupation ( a non- 
zero density in the thermodynamic limit) [26]. Introduced in [!)] in the context 
of an ideal Bose gas, it was due to naturally arising interactions a difficult 
problem to realize Bose-Einstein condensation experimentally [8, 20]. 
As we will see, also the formal description is more elaborated. To motivate the 
definition of BEC for vanishing temperature we follow the continuous approach 
in [26]. To formalize the concept of a macroscopic occupation of a single particle 
state we remember the definition of the one-particle density matrix [26], i.e. the 
operator on L^(M^) given by the kernel 

N 

7(x, x') = N X2, . . . , xn)^^^{x', X2, . . . , xn) n ^^i 

with the normalized ground state wave function of the many Boson Hamil- 
tonian. BEC in the ground state is then defined that the projection operator 7 
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has an eigenvalue of order N in the thermodynamic Hmit. 

Remembering that for the ideal Bose gas the multi-particle ground state can 
be represented as a product 



of the single particle ground state ipo the one-particle density matrix becomes 

thus the definition of BEC above is natural and can also be related to the 
thermodynamic formalism (see e.g. [24, 26] and references). In particular, it is 
of interest to consider BEC for the ideal Bose gas with a random background 
potential. In this case the Lifshitz tail behavior at the bottom of the spectrum 
makes a generalized form of Bose-Einstein condensation possible even for d = 1, 
2 (see [24] and references cited there). 

The situation in the Gross-Pitaevskii-limit is close to the situation for the ideal 
Bose gas [26]. The one-particle density matrix is asymptotically given by 



Physically the content of (7) is that all Bose particles will condensate in the 
GP ground state motivating the definition of complete (or 100%) BEC in [26]. 

The purpose of the present publication is a first step to analyze the fine struc- 
ture of the Gross-Pitaevskii ground state. Under the assumption of a random 
background potential we want to understand how tp'^^ is related to the eigen- 
states of the single particle Hamiltonian. More familiar is this problem in the 
following two settings. 

If the Bosons are trapped by a potential tending to oo, i.e. liminf V{x) = 
oo, the spectral properties of the single particle are invariant in the thermo- 
dynamic limit, i.e. the discrete spectrum and the strictly positive distance 
between the first two eigenvalues. Assuming Na — t- in the continuous setting, 
respectively NU in the context of the discrete Gross-Pitaevskii model, the 
interaction energy is a small perturbation of the single particle energy func- 
tional. In this situation it is natural, that in the thermodynamic limit ip*^^ and 
the single particle ground state ipo coincide [25]. 

A complementary situation is given if the Bosons are confined to a cube A^ with 
I Ail — )• oo but without a background potential. As described in [26] assuming 
p = N/ LP' and g = Na /L in the limit N ^ oo one can prove 



$Q^^(xi,...,X^)=a=l MX^) 



(7) 



7(x,x') 



N-^oo 



iV(^GP(x)<^GP(x'). 




7(2;, y)dxdy = 1 



i.e. BEC in the normalized single particle ground state po = L ^^^^XA^- As 
explained in [26] 



Na _ pa 



L 1/L2 
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is in this context the natural interaction parameter since " in the GP hmit the 
interaction energy per particle is of the same order of magnitude as the energy 
gap in the box, so that the interaction is still clearly visible". 

As emphasized in the physics literature (see e.g. [4, 27]), new phenomena like 
fragmented BEC (Lifshitz glasses) should occur when Bosons are trapped in a 
random background potential. Our purpose in this publication is more modest. 
We want to understand the natural interaction parameter in a random media, 
s.t. the Gross-Pitaevskii ground state is close to the ground state of the single 
particle Hamiltonian as it is suggested by the situation in the ideal Bose gas. 
As we will see the setting of Bosons trapped in a random potential is not really 
comparable to the two situations described above. 

Under our assumptions, near which is almost sure limit of mf{H^ we are 
in the localized regime, i.e. one has pure point spectrum and localized eigen- 
functions. In contrast to the situation with vanishing potential the eigenstates 
close to the bottom of the spectrum are localized in a small part of A^, i.e. 
the interaction energy will be larger than in the case of the homogeneous Bose 
gas. In the random case, we determine the almost sure behavior of the ground 
state from information on the integrated density of states (see Lemma 5). Un- 
der our weak Lifshitz tails assumption (H3), we obtain that the ground state 
energy is of size (logL)"^/*^. When L — >■ +oo, the difference between the first 
two eigenvalues will tend to zero; the speed at which this happens is crucial 
in our analysis (see Proposition 10). In our case, we estimate that, with good 
probability, it must be at least of order L~'^. This difference is much smaller 
than the one obtained in the homogeneous Bose gas where it typically is of or- 
der L~^. We deem that the estimate L~'^ for the spacing is not optimal in the 
present setting. This estimate is the correct one in the bulk of the spectrum; 
at the edges, the spacings should be larger. It seems that getting an optimal 
estimate requires a much better knowledge of the integrated density of states 
or, in other words, much sharper Lifshitz tails type estimates (see (113)) and 
Minami type estimates that take into account the fact that we work at the edge 
of the spectrum (see (112)). Combining these observations explains the interac- 
tion parameter U = o(L~'^/i^^ (log L)) that we don't believe to be optimal. 

Let us now briefly outline the structure of our paper. To prove our result we 
need two ingredients. We need an upper bound of the interaction term, i.e. we 
have to estimate the ||.||4- norm of the single particle ground state (po. At the 
same time, we need a lower bound of the distance of the first two single particle 
eigenvalues asymptotically almost surely (a.a.s.) i.e. with a probability tend- 
ing to 1 in the thermodynamic limit. Comparing these two estimates we will 
see that under the assumptions of Theorem 3 it is energetically favorable, that 
the Gross-Pitaevskii ground state and the single particle ground state coincide. 
This will be proven at the end of this publication. 

To estimate the interaction term we will prove in Lemma 5 that almost surely 
in the thermodynamic limit the single particle ground state is fiat, i.e. 

ii_ ||2 L— ^oo 

llVy^oll > U a.a.s. 
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This then yields an estimate of the interaction term which is the purpose of 
Proposition 4. 

The a.a.s. lower bound of the distance of the first two single particle eigenvalues 
is a little bit more intricate and uses the Wegner and Minami estimates; it is 
related to the methods developed in [: 2]. In Lemma 12, we first estimate 
the probability that the first two eigenstates and also their localization center 
are close together. If the localization centers are relatively far away, one can 
decouple the eigenstates and treat the first two eigenvalues of each other. This 
is used in Lemma 13. 

2. Estimating the interaction term 

The main result of this section is an upper bound on E^'^ - E^[a;,L]. This 
quantity is non negative (see (2)) and we prove 

Proposition 4. There exists C > 0, such that, for any p S N, one has 

(8) P[eSP -EP[^,L] <C[//,(logL)] >l-L-f 
where fd is defined in (3). 

By definition, for ipo{uj, L) the ground state of one has 

(9) E^P < e^M^^L)] = EP[a;,L] + U\\^^{u:,L)\\\ 
resp. 

E^P -EPKL]<[/||(^oK^)||!. 

To prove Proposition 4, resp. control the interaction term, we first estimate the 
ground state energy of the random Schrodinger operator and derive in Corollary 
6 an estimate on the "flatness" of its ground state. We start with the Dirichlet 
and Neumann boundary cases. 

Lemma 5. Assume (H3) is satisfied. Let Eq{uj,L) be the ground state energy 
of and denote by ipo{uj,L) the associated positive normalized ground state. 
Then, for any p > 0, there is a constant C > such that, for L sufficiently 
large, 



(10) P \c-\logL)-^/'^ < E^{uj,L) < C{\ogL)- 



2/d 



> 1 - L-P. 



As K; is non negative and (po{(^, L) normalized, one has \\'S/lpq{uj, L) |p < Eq{uj, L). 
Hence, Proposition 4 implies the following "flatness" estimate of the ground 
state. 

Corollary 6. Under the assumptions of Proposition 4, for any p > 0, there is 
a constant C > such that, for L sufficiently large. 



(11) 



||V<^o(^,i)f <C(logL)-2/'^' 



> 1 - L-P. 
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It is maybe interesting to note that from a Lifshitz tail type estimate (i.e. the 
annealed estimate), we recover the (approximate) almost sure behavior of the 
ground state energy of ^ (i.e. the quenched estimate) (see e.g. [35]). 
We note that Proposition 4 and Corollary 6 also hold if we replace the periodic 
ground state and ground state energy by the Neumann or Dirichlet ones. 

Proof of Lemma 5. Fix i>l. Decompose the interval [—L,L] into intervals of 
length comprised between i/2 and 2i. This yields a partition of in paral- 
lelepipeds i.e. 

Al= U 

such that 

• Pj = Ij X ■■■ X If where the intervals {Ij)i<k<d satisfy £/2 < < 2i 
. for j / /, Pj n Pj, = 0, 

• J, the number of parallelepiped, satisfies 2-'^{L/lY < J < 2'^{L/ef. 

In the continuous model, one can take the parallel piped to be cubes. 

Denote by oo\\^ the restriction of a; to A^. Furthermore, let lo^'^ be the periodic 

extension of i^ia^ to 'E'^ i.e. for /? € A^^ and 7 € Z"^, wj^— = ujjj where 

L = 2L + 1 in the discrete case and 2L in the continuous one. As is 
the periodic restriction of i?^ to Ai, we know that EQ[a;,L] = inf a{H^p.L) 
where this last operator is considered as acting on the full space M'^ or Z'^ (see 
e.g. [19]). 

We can now decompose M'' or Z'^ into U^g^d U^=i {'jL + Pj). By Dirichlet- 
Neumann bracketing (see e.g. [15, 17]) H^p,l satisfies as an operator on Z'^ or 

(12) ®-yezd ®'j=i^u\(^L+p,) - ^^^'"^ - ®/=i ^uli-rL+Pj)- 
Define 

(13) Ej[w,£,j] =infa(H'|pP for •e{N,D}; 

here, the superscripts D and N refer respectively to the Dirichlet and Neu- 
mann boundary conditions. As uj^'^ is LZ'^-periodic, H'ln and H*,, -p „ , are 

unitarily equivalent. The bracketing (12) then yields 

inf EN[a;,£,j] <EP[u;,L] < inf E°[u;,^,j]. 

Labeling every second interval of the partition of [—L, L] used to construct 
the partition of A/,, we can partition the interval {1, • • • ,J} into 2^^ sets, say 
(J7i)i<;<2ti such that 

(1) if / / Ji n Ji, = 0, 

(2) for j e Ji and / € Ji such that j / j', one has dist(Pj, P^/) > i/2, 

(3) there exists C > such that for 1 < / < 2'^, C-^{L/lf < #Ji < 
C{L/iY. 
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Assume R is given by (HO). By (2) of the definition of the partition above, for 
any / > 2R, all the {Rl^p^)j^j^, resp. all the {E'o[u;,i,j])j^j^ (for • G {N,D}) 
are independent. Hence, using (13), we compute 

P[EP[a;,L] > E] < P[inf eD[^, ^, j] >E]<Y,ll F[E?[u;,^,j] > E] 
^ 1=1 j&Ji 

=En(i-'p[E?k^,j]<i?]). 

Pick E = ci~'^ where c is given by assumption (H3) and 

{k log L - c^^ log c)^/'^ < ^ < (A: log L - c"^ log c)^/*^ + 1 
where k will be chosen below. Applying the Lifshitz estimate (H3), we obtain 

2"^ 

P[EP[a;,L] > i?] < J] (l - e-^i°s ^^^^ 
1=1 

2d 

z=i 

if we choose k < cd as C'^L/if < #Ji < C{L/ef for I < I < 2'^. 
Hence, we have 

P[E^[w,L] <E] > l-0(L-°°). 
To estimate from below, we use again (13) to get 

F[E^[uj,L] < E] < J]P[E^[u;,^,j] <i?]). 

Pick E = Ci~'^ where C is given by assumption (H3) and 

{k log L - log C)^/'^ < ^ < (A: log L - C'^ log C7)i/'^ + 1 

where k will be chosen below. As ^J' < C{L/lY, applying the Lifshitz estimate 
(H3), we obtain 

if we choose k > {d + p)C. Hence, we have 

P[E^[a;,L] > ^] > l-L'P. 
This completes the proof of Lemma 5. □ 

To prove estimate (8), we will use the spectral decomposition of — A£. Though 
the arguments in the discrete and continuous cases are quite similar, it simplifies 
the discussion to distinguish between the discrete and the continuous case rather 
than to introduce uniform notations. We start with the discrete case. 
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Lemma 7. There exists C > such that, for e G (0, 1) and L € N satisfying 
L-e>l one has for u G f{'L'^/{2L + 1)1'^) with \\u\\2 = 1 and (-Afn, u) < 
the estimate 

i/(i<3, 

(14) ||u|[4 < Cgd,{e) where gd{() = I ^llog^l if d = A, 

_e ifd>5. 

Proof. The spectral decomposition of — A£ is given by the discrete Fourier 
transform that we recall now. Identify with the Abelian group Z'^/(2L + 
For u G Ul, set 



(15) u = iu^)\^\<L where = ^ ^^^^^ ^ 
Then, one checks that (see e.g. [21]) 



up ■ e 



-2i7r7/3/(2L+l) 



m<L 



(16) (— A£m)" = (/i(7)ii7)|'y|<L where h{'y) = 2d — 2 cos 



2L + 1 



Pick u € r(Z''/(2L + 1)Z'*) with ||n||2 = 1 and (-A£u,ti) < and write 
u = J2k=o ^fc where fc^ G N, — log e < /cg < — log e + 1 and 

• Uq = U- l|^|<eL 

• for 1 < /C < /Ce - 1, % = U • le'=~ieL<|7|<e'=£L 

where u denotes the discrete Fourier transform defined defined in (15). 
Then, for k ^ k' , {uk,Uki) = and, using (16), for /c > 1, 



Uk,Uk) 



-Afn, n) < 



A;=0 



k=0 



I.e. 



(17) 



k=0 



Hence, using (15) and Holder's inequality, we compute 

. . 1 



KUkjl3\ 



< 



< 



{2L + 1)^^/2 
1 

(2L + l)'^/2 
1 



e'=-ieL<|7|<e*eL 
i^e'=-ieL<|7|<e'=eL 

\ukh {eh{2L + l)Y/^. 



2i7r7/3/(2L+l) 
1/p 



l/<? 



E 1 

^e''-ieL<|7|<e'''eL 



(2L + l)'^/2 
So, for p = (7 = 2, one gets 

(18) \\uk\\oo<\\ukh{e^eY'\ 
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Then, using (17), we compute 

k k k 

\\u\U < J2 ll^'^ll^ ^ E VhkhWukWoo < cX;e'=('^-4)/4e'^/4 < Cg^ie) 

k=0 k=0 k=0 

where gd is defined in (14). This completes the proof of Lemma 7. □ 

Remark 8. Lemma 7 is essentiaUy optimal as, for L sufficiently large, 

• the trial function 

^ feif7 = 0, 

\(2L + l)-'^/2if7/0, 

satisfies 1 < ||n||2 < 1 + e, (— A^n, u) < Ce^ and ||n||4 > e/C; 

• the trial function 

f (2eL + l)-'^/^ if |7| < eL, 
^ [0 if l7| > eL, 

satisfies ||n||2 = 1, {—Af^u,u) < Ce^ and |ln||4 > e'^/^/C 
We now turn to the continuous case. 

Lemma 9. Fix ij G (0,1/4). There exists C > such that, for e G (0,1), 
n > {d-2)ri~^ + l andL^n satisfying L-e > 1 one has for u G H"'{R'^/{2L)Z'^) 
with {—Af^u,u) < the norm estimate 

fC'^/^ ifd<3, 
C^-^llogel ifd = A, 



Wh < C gd,n{£)\\u\\'Hn where gd,r){0 



ifd>5. 



Proof. We now use the Fourier series transform to decompose — A|^. Identify 
Al with the Abehan group R'^/2LZ'^. For uGTLl, set 

(19) u = (n^)^g^. where u, = u{e) ■ e~^^^'/^de. 
Then, 

(20) <e) = ' -7^-^^^/" 



and 



(21) {-A^^uY 





TT-J 




L 



if u G Vl. 

|7|<L 



Pick u as in Lemma 9 and decompose it as in the proof of Lemma 7 i.e. write 
u = Ylk=o where A;^ G N, — log e < < — log e + 1 and 

• Uq = U- l|7|<eL 

• for 1 < k < ke - 1, Uk = U ■ le'=-leL<|7|<e'=£L 

• Uk^=U - lefcEeL<|7| 
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where u denotes the Fourier series transform defined in (19) and (20). 

The control on Uk for < A; < fc^ — 1 is obtained in the same way as in 

the proof of Lemma 7 namely the estimate (18) holds for < A; < — 1. The 

additional ingredient that we need is to obtain a control over the large frequency 

components. 

Recall that 



l7/l|2 



L 



n/2 



Fix r > d. For notational convenience, write v = and compute 



1 



< 



(2L)'^/2 
1 

(2Lp72 



E ( 

e^i^ eL<.\'y\ 

E 

e'^s eL<\'y\ 



7r7 



L 



r/2 



\(yy 



7r7 



-r/2 



< 



(2L)'^/2 



E 



e*-'~l£L<|7|<e*-'£L 



7r7 



L 



^e'=-ieL<|7|<e'=eL 

<C\\v\\l^-e^-^ 




r-2/q 



i/(2g) 



1/2 



li p = T], q = 1 — Tj and r = (n — 1)77 + 2>(iasn> (d — 2)77 ^ + 1. One then 
completes the proof of Lemma 9 in the same way as that of Lemma 7. □ 



Proof of Proposition 4 ■ In the discrete case Proposition 4 is a consequence of 
Lemma 5 and Lemma 7 with e = C {log L)~^ / . 

To be able to apply Lemma 9 to (po{uj,L) in the continuous case, we need to 
show that, for any n > d, ipo{u},L) E with a bounded depending only on 
n not on L or w. Therefore we use the first assumption on the random field 
Vuj i.e. that, for any a € N'^, ||9"K;||a;,a;,oo < +00. Hence, as ip(){uj,L) is an 
eigenvector of — A£ + V^, using the eigenvalue equation 

inductively, we see that 

||||</'o('^,-^^)||j/"L,oo < +00. 

Proposition 4 in the continuous case is then a consequence of Lemma 5 and 
Lemma 9 with e = C(log L)~^^^. This completes the proof of Proposition 4. □ 
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3. The spectral gap of the random Hamiltonian 
The main result of the present section is 

Proposition 10. Let the first two eigenvalues ofa^^ be denoted 62/Eg [a;,L] < 
E^[a;,L]. Then, forp > 0, there exists C > such that, for L sufficiently large 
and rj € (0, 1), one has 



Ef[uj,L] - E^[oj,L] < r]L^'^] <C7]\l + (log L)'^-2/d+e 



with e = in the discrete setting resp. e > arbitrary in the continuous case. 

In the localization regime, both the level-spacing and the localization centers 
spacing have been studied in e.g. [12, 14]. The main difficulty arising in the 
present setting is that the interval over which we need to control the spacing 
is of length C(log L)"^/*^; it is large compared to the length scales dealt with 
in [12, 14]. 

Our analysis of the spectral gap relies on the description of the ground state 
resulting from the analysis of the Anderson model i/^ in the localized regime 
(see e.g. [15], [33]). Under the assumptions made above on H^^, there exists / a 
compact interval containing such that, in I, the assumptions of the Aizenman- 
Molchanov technique (see e.g. [1, 2]) or of the multi-scale analysis (see e.g. [11]) 
are satisfied. One proves 

Lemma 11 ([11, 22]). There exists a > such that, for any p > 0, there exists 
q > such that, for any L > 1 and ^ G (0, 1), there exists ^i^s,L C 0, such that 

. P[0/,5,l] > 1 - L-P, 
• foruj G J^/,<5,L! 

one has that, if<fin,Lu 'is 0, normalized eigenvector o/H(^|y^^ 
associated to En^uj £ I, CLnd Xn{oj) € is a maximum of x ^ |(/5„^^(a;)| 
on Kl then, for x G A^, one has, 

I , ^„ I e~"l^~^"('^)l in the discrete case, 

(22) \ipn,Ux>\ < ■ < „„u_^ (^)|« . 

I g .Lnyujjf continuous case. 

Note that, for a given eigenfunction, the maximum of its modulus need not 
be unique but two maxima can not be further apart from each other than a 
distance of order logL. So for each eigenfunction, we can choose a maximum 
of its modulus that we dub center of localization for this eigenfunction. 
To prove Proposition 10, we will distinguish two cases whether the localization 
centers associated to Eq[w,L] and Ej'[ti;,L], say, respectively xo{uj) and xi{lo) 
are close to or far away from each other. 

In Lemma 12, we show that the centers of localization being close is a very rare 
event as a consequence of the Minami estimate. 

In Lemma 13, we estimate the probability of Eo[t(;,L] and Ei[ct;,L] being close 
to each other when xo{uj) and xi{uj) are far away from each other. In this case, 
Eo[a;, L] and Ei [w, L] are essentially independent of each other, and the estimate 
is obtained using Wegner's estimate. 
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Lemma 12. For p > 0, there exists Lq > such that, for X > 0, L > Lq and 
rj G (0, 1), one has 

"En^,L] -EP[l^,L] <r?L-^ 

_|xo(a;)-xiH| < A(logL)i/«_ 

with ^ = 1 in the discrete setting resp. ^ > 1 arbitrary in the continuous case. 



Proof. Let us start with the discrete setting. Fix p > and let q be given 
by Lemma 11. The basic observation following from Lemma 11 is that, for 
uj G J^7,5,L) if Xnii^) is the localization center of (fni^, L) and I < L, then 

(23) ||(/?°-E^(u;,L))^„(a;,L,/)|| + |||(^„(a;,L,/)|| -1| <CL«e-'. 



where 



= [H^ L\\xnH+Ai is H^L restricted to the cube Xn{uj) + A/, 



• ipn{^,L,l) = lx„(a;)+A;<^n('^, is the eigcnfunction ipn{uj,L) restricted 
to the cube Xn(u;) + A;. 

To apply the observation above we pick a covering (C'j)o<j<J of A/, by cubes 
of side length of order logL i.e. Al C Uo<j<j^j- Then the number of cubes 
J can be estimated by J < CL'^(log and there exists C > (depending 
on A, q and v) such that, if \xq{uj) — xi(ti;)| < AlogL and / > CAlogL, there 
exists a cube Cj (containing xq{uj)) such that, for L sufficiently large 

1 

J2 {mi - E^{u,L))^k{^,L,j)\\ + |||(^,(a;,L,i)|| - 1|) 
fe=o 

+ \{M^,L,j),M^,LJ))\ <L"72 

where we have set q — CXa < —u (see (23)) and 

• HIj is the operator restricted to the cube Cj + A;, 

• 'fk{^,L,j) = lcj+Ai(pk{^,L) for k G {0,1}. 

Let C be given by Lemma 5 and define / = [0,2C{log L)~'^^^]. Decompose 
/ C U^i+^/^ where 

• Im are intervals of length 4r]L~'^, 

• for m G {0, . . . , M - 1}, hm H /2(m+i) = = hm+i n hm+S, 

• for m G {0, . . . , M}, n hm+i is of length 2r]L~^. 

One can choose M < CL'^(log L)~^/'^7y~^. This implies that, for L sufficiently 
large, 

r Ena;,L]-EP[a;,L]<r?L-'^\ 
\xo{u}) - xi{uj)\ < AlogL 



where = \ ^li^s,L and 



J 2M+1 

^2 = U U {{Huj)\Cj+Ai has two eigenvalues in 

j=l m=0 
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By Lemma 11, we know that 

P[f]i] < L-'P 

Minami's estimate (H.2) and the esthnate on M tells us that 

P[f]2] < C L^'^ {\og L)-'^''^''^ri-^{i]L-'^{C log lYf < Cr]{logL f-^/'^. 

This completes the proof for the discrete setting. The proof for the continuous 
case is very similar. One has to replace lx„(aj)+Aj by a smooth version of the 
characteristic function of the cube Xn(c<j)+A; (see for example [34]), resp. change 
the length scale logL to (logL)^/'> in the side length of the boxes where one 
restricts the eigenfunctions. This is necessary because of the weaker estimate 
in Lemma 11. This completes the proof of Lemma 12. □ 

We now estimate the probability of the spectral gap being small conditioned on 
the fact that the localization centers are far away from one another. We prove 

Lemma 13. For any p > 0, there exists A > and C > such that, for L 
sufficiently large and rj € (0, 1), one has 

■Ef[w,L] -E^[a;,L] < r^L^^, ' 
Jxo(a;)-xi(a;)| > A(logL)V«_ 

with ^ = 1 in the discrete setting, resp. > 1 arbitrary in the continuous case. 



<Ci] + L- 



Proof. Using the same line of reasoning as in the proof of Lemma 12 we give 
the proof in the discrete setting. 

Fix u > 2d + p and split the interval [0, C(log L)~^/'^] into intervals of length 
L~'^ as in the proof of Lemma 12. By Minami's estimate, we know that, for L 
sufficiently large 

(24) F[E^[uj,L]-E^[u},L]<L-''] < C(log L)-2/'^L^l2('^"^) < . 
So we may assume that E5^[ct;,L] — Eq [a;,L] > L~'^. 

As in the proof of Lemma 12, pick a covering of Al by cubes, say {Cj)o<j<j 
of side length less than logL such that J, the number of cubes, satisfies J < 
C{L/\ogLf. 

Assume that Cj is the cube containing 2:1(0;), Ei[a;,L] — Eo[w,L] < i]L~'^ and 
\xq{uj) — xi{oj)\ > AlogL. Let Aj = Aj;, \ (Cj +A3/4;^ log l). Define the operators 
{H^)\l^c, resp. {Huj)\Cj+k^i„^L/A *° restriction of i?^^ to A^, resp. Cj + 

■^AiogL/4) with Dirichlet boundary conditions. If A > 8 and L is large enough, 
we know that 

• dist(xo(t^),5Ap > AlogL - 3/4AlogL - logL > AlogL/8, 

• disi{xi{uj),d{Cj + Ax\o^L/i)) > AlogL/4 
. dist(A^-,Cj + A;,iogi/4) > AlogL/2 > R 

with > as in the decorrelation assumption (HO). Hence, for A sufficiently 
large, using the estimate (23) for the operators {Hi^)\f^c and 
we know that: 
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• The operator {Hic)\Cj+KxiogL/i a-dmits an eigenvalue, say Ei{uj), that 
satisfies \Ei{uj) - E[{uj)\ < L"^"; 

• The operator (ffi^)|^c admits an eigenvalue, say Eq{uj), that satisfies 

j£'o(a;) — Eq{uj)\ < L~'^'^ . Moreover, as {H^)\jyc^ is the Dirichlet restric- 
tion of ^, its eigenvalues are larger than those of ^. In particular, 
its second eigenvalue is larger than £^1(0;). Hence, up to a small loss in 
probability, we may assume it is larger than Eq{uj) + L~'^ as we know 
the estimate (24). This implies that we may assume that Eq(u}) is the 
ground state of {H^^)\c. 

So we obtain 

;^Pr, . Ti T^Pr, . Ti / „ T-d ^ 



to: 



E{[uj,L] - E^[u},L] <r]L- 



1<3<J 



|xo(f^) — 2;i(a;)| > AlogL 
where VLi=VL \ U {w; Ei[a;,L] — Eo[a;,L] < L~'^} and 

% = {^;dist(cT((i7^)|c,+A,,g,/J,infa((F^)|Ac)) < ^L"" + L~^] 

As {Hu))\Cj+Kx\osL/4 ^^"^ {Hu])\K'=. are independent of each other, we estimate 
the probability of using Wegner's estimate to obtain 

P[%] < C{r]L-'^ + L-''){\ogLY. 

Hence, one obtains 

^ ' ^ 01 , J _ / , < c(^-'^ + 2.-n(logL)'^- -T + 2L-P 

|xoM-xiH| > AlogL J" ^ ^ (logL)rf^ 

< C(r/ + L-P) 

\i V > d + p. 

This completes the proof in the discrete setting. To prove Lemma 13 for the 
continuous case, one does the same modifications as in the proof of Lemma 12 
in the continuous setting. □ 

Setting e = d{l/^ — 1), Proposition 10 then follows from Lemma 12 and 
Lemma 13. 

4. Proof of Theorem 3 

Defining ttq = |99o)('/'o| and applying the definition of the ground state, we can 
estimate 

L] II (1 - ^o)</'^^ II + EP[a;, L] |KoV^^^ II 

^T?GP||/i ^ \, GPii , -pGP 11^ ,„GP|| 

respectively 

(Ena;,L] - E^P) ||(1 - vro)^^^! < [^Z " ^l[^M) Iko^^^^H- 
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As a consequence of Proposition 4 and Proposition 10, we know with a proba- 
bility larger than 1 — (Cr/ + L~p) that, for r] G (0, 1) and fd defined in (3) the 
estimates 

Enu;,L] - > Ef[u;,L] - eJ- [a;,L] > rjL-''[l + (log L)^-2/'^+r ' 

and 

E^P -EP[a;,L] <CC//,(logL) 

are satisfied. We obtain 

11(1 - 7ro)(^«P|| < CC//rf(logL)7?-^L'^[l + (logL)'^-2/'^+^] ||7ro(^GP|| 

and 

|(<,o,^^^)P = ||vroc^^^f = l-||(l-vroV^^f 

> 1 - [(:7[//d(logL)r/-iL'^[l + (logL)'^-2/d+e] 
Applying the assumption concerning the coupling constant U i.e. 

U = U{L) = o (l-'^[1 + {logLf-'/''+r'mogL)]-'^ 
and setting 

r] = riiL) = ^\U{L)Ld\l + (logL)'i-2/d+6]/,(iogL)| 
we get that, ri{L) — )■ when L +oo and for some C > 0, 

P({a;; ||((^o, V^^^)! - 1| > Cii{L)}) < C{v{L) + L'^). 
This completes the proof of Theorem 3. 
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